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Abstract—  Scattering  from  the  seafloor  Is  often  assumed  to  be 
controlled  by  two  spatial  scales:  the  larger  scale  associated  with 
reflections  from  plane  facets  and  the  smaller  one  associated 
with  diffuse  scattering  from  height  variations.  Choosing  the 
wavenumber  for  this  partitioning  has  proven  to  be  Important  but 
troublesome.  For  this  work  scattering  data  Is  simulated  using 
Helmholtz-KIrchhoff  or  physical  optics  theory  and  selected  Input 
geomorphology.  These  data  are  Inverted  using  a  simpler  two- 
scale  roughness  model.  The  work  Introduces  a  new  criterion 
for  effecting  this  partition  based  on  setting  a  roughness 
parameter  equal  to  unity.  The  criterion  Is  shown  to  be  valid  for 
the  cases  analyzed  based  on  the  ability  of  the  Inversion  model 
to  recover  the  Input  geomorphology. 


I.  INTRODUCTION 

Among  the  approaches  to  the  study  of  rough  surface 
scattering  two-scale  roughness  (TSR)  theory,  also  called 
composite-roughness  theory,  is  particularly  attractive 
because  it  separates  scattering  into  what  appears  to  be  its 
two  primary  mechanisms:  scattering  controlled  by  large 
random  facets  and  scattering  controlled  by  small 
roughness  1 1,2].  The  reported  advantage  of  this  approach 
is  that  its  applicability  can  span  a  wide  range  of  surface 
roughness,  since  it  allows  for  the  selective  application  of 
appropriate  approximations  in  scattering  theory.  While 
this  is  certainly  true,  it  leads  one  to  attempt  to  answer 
the  troublesome  question  of  where  to  set  the  scale 
partition  based  on  the  requirement  for  the  validity  of  the 
approximations.  The  question  should  rather  be  posed  in 
terms  of  "At  which  scales  is  the  scattering  like  reflection 
from  a  collection  of  facets  and  at  which  scales  is 
scattering  like  diffraction  from  roughness?”.  The  goal  of 
this  work  is  to  allow  the  scattering  phenomenon  itself  to 
suggest  what  should  be  the  scale  partition. 


This  work  is  supported  by  the  Acoustic  Reverberation  Special 
Research  Program  (ARSRI*)  of  the  Office  of  Naval  Research  (ONRy 


II.  APPLICABLE  SCATTERING  MODELS 

Here  we  introduce  the  term  Rough  Facet  Model 
(RFM)  to  refer  to  a  part  of  a  broader  set  of  programs 
called  the  Bistatic  Scattering  Strength  Model  (BISSM) 
[3].  RFM  is  that  part  that  is  the  mathematical 
description  of  the  scattering  mechanism  that  includes  a 
product  between  a  factor,  F(» ),  describing  fine-scale  slope 
scatter  and  a  factor,  M(*),  for  scatter  from 
microroughness.  RFM  allows  inversion  of  the  scattering 
problem  and.  thereby,  the  extraction  of  slope  and  height 
roughness  for  the  scattering  surfaces.  Since  no  reliable 
experimental  data  are  available,  scattering  data  arc 
simulated  using  a  high-fidelity,  2-D  scattering  model  and 
selected  input  seafloor  morphology.  This  high-fidelity 
model  is  based  on  Monte  Carlo  of  evaluation  of  the 
Helmhollz-Kirchhoff  (H/K)  integral  over  an  ensemble  of 
surfaces  with  the  desired  statistical  properties.  RFM 
predictions  are  fitted  to  results  from  the  H/K  simulations 
and  the  scale  partition  is  then  established  when  successful 
recovery  the  correct  geomorphology  is  attained. 

The  RFM  algorithm  gives  backscatter  at  near  normal 
incidence  as 

mbs  =  M(k,V9).  F(5f,0)  (1) 

where  k  is  the  acoustic  wavenumber;  is  the  standard 
deviation  of  the  microroughness;  6f  is  the  standard 
deviation  of  the  fine-scale  slopes;  and  8  is  the  angle  of 
incidence  measured  from  normal  (zenith  angle).  For 
convenient  visualization  we  define  and  refer  to  the  slope 
angle  af  =  tan‘!6f. 

The  coherent  reflection  coefficient  is  given  by 

M  =  cxp[-gj  (2) 

where  g  is  the  well  knowyi  roughness  parameter  restricted 
here  it  to  the  microroughncss,  i.c.,  g^  =  k*cos^8- 
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(5) 


The  function  F(»)  is  the  high-frequency  limit  of  the  H/K 
theory  in  the  form  derived  Brekhovskikh  and  Lysanov  [4]: 


F  = 


R 


2 

o 


. exp 

Sax  6  f'cos^B 


r  tan^0 


(3) 


where  R0  is  the  Rayleigh  reflection  coefficient.  For  this 
work  R„  =  1. 


6  f2  =  2tt 

In  Eqs.  (4)  and  (5),  K(  and  Kh  indicate  the  low-  and 
high-frequency  wavenumbers  present  on  the  surface, 
respectively.  These  wavenumbers  are  determined  by  the 
size  of  the  footprint  and  the  grid  spacing,  respectively. 
Kc  denotes  the  partition  (cut-off)  wavenumber. 


Kc 

A2(K)K3dK 

K, 


To  simulate  scattering  data  we  use  a  model  that 
consists  of  evaluating  numerically  the  full 
Helmohltz-Kirchhoff  integral  with  a  minimum  number  of 
approximations.  A  complete  description  of  the  technique 
can  be  found  in  ref.  5. 


III.  NUMERICAL  SIMULATIONS 


Although  this  is  a  purely  numerical  study,  a  practical 
scenario  based  on  a  multibeam,  swath  bathymetric  sonar 
system  is  adopted  here  because  it  will  allow  easy  testing 
of  the  results  and  because  we  envision  that  system  as  a 
major  application  of  these  results  once  they  are  validated. 
However,  the  scenario  is  also  applicable  to  the  low- 
frequency  active  sonar  system  used  in  the  Acoustic 
Reverberation  Special  Research  Program  (ARSRP)  of 
the  Office  of  Naval  Research  (ONR).  We  apply  it  here 
at  frequencies  more  relevant  to  ARSRP  data-frequencies 
equal  to  250,  500,  1000,  and  2000  Hz.  The  water  depth 
is  chosen  to  be  2000  m  and  the  beamwidth  is  two 
degrees.  Zenith  angles  are  set  at  3  to  43  degrees  in  two- 
degree  increments. 


Ensembles  of  2-D  (i.e.,  surface  heights  are  function  of 
x  and  y)  randomly  rough  surfaces  were  generated  using 
a  filtering  technique  developed  by  Caruthers  and 
Novarini  [6,7],  Amplitude  spectra  used  to  describe 
seafloor  roughness  is  of  the  isotropic  power-law  form, 
A(K)  =  aK'b,  where  K  is  the  wavenumber  of  the  surface 
roughness  (in  cycles/meter). 


It  should  be  noted  that,  in  the  RFM  scattering  model, 
the  two  surface  parameters  (the  rms  height  of  the 
micro-scale  surface  and  the  rms  slope  of  the  fine-scale 
surface)  arc  band  limited  quantities.  The  variance  of  the 
microroughncss  for  a  2-D  isotropic  surfaces  is  given  by 


v  =  2* 


K„ 


A^(K)KdK 


(4) 


K„ 


The  variance  of  the  slopes  of  the  fine-scale  surface  by 


Table  I  shows  the  spectral  parameters  a  and  b  for  the 
surfaces  used  in  this  work.  Given  in  the  table  areoT,  the 
total  rms  height  (i.e.,  presuming  that  all  the  deformation 
of  the  surface  is  considered  to  be  microroughness)  and 
ST,  the  total  rms  slope  (i.e.,  presuming  that  all  the 
deformation  of  the  surface  in  attributed  to  slope 
variance).  Also  given  in  the  table  is  the  slope  angle  aT. 
The  dimensions  for  each  of  the  surfaces  is  256  by  256 
meters  (1024  x  1024  grid  points,  with  a  grid  spacing  of 
0.25  m).  Each  surface  realization  was  split  into  a 
number  of  ensonified  regions  for  ensemble  averaging 
purposes.  The  number  of  subregions  depends  on 
geometry  of  the  simulation  and  the  zenith  angle  of 
incidence. 


TABLE  I.  PA  RAMETERS  OF  'll  tE  POWER-LAW  SPECTRA  FOR 
THE  SIMULATED  SURFACES 


surface 

label 

spectra!  parameters 
a  b 

OF 

2.32xl02 

1.75 

01 

1.86x1a2 

1.75 

0D 

6.00xia3 

2.27 

01 

4.80xI03 

2.27 

0J 

3.o  xia3 

2.27 

rms  values 


<Mm) 

«T 

<*T  (*8) 

1.80 

0.594 

30.7 

1.45 

0.477 

25.5 

4.47 

0.464 

24.9 

3.58 

0.372 

20.4 

224 

0.233 

13.1 

IV.  RESULTS 

Here  we  hypothesize  that  the  scattering  phenomenon 
itself  will  establish  a  facet  size  such  that  the 
microroughness  that  exists  upon  it  leads  to  g  »1. 
Larger  facets  contain  too  much  roughness  to  produce  a 
coherent  reflection.  Smaller  facets  arc  not  the  full  extent 
of  the  region  that  can  cause  a  coherent  reflection. 
Therefore,  the  partition  wavenumber  can  be  determined 
from  Eqs.  (4)  and  Sf  determined  from  Eq.  (5).  For 
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simplicity  we  will  approximate  gM  for  all  angles  (0  to  45 
deg)  by  its  value  at  9  =Cf ,  i.e.,  g/10  =  4cr#42k^s  1. 


Fig.  1.  Comparison  of  the  RKM  and  H/K  results  at  500  Hz  for 
surface  0D  for  different  values  of  spectral  slope  b  assuming  that 

W  =  10' 


Adopting  the  forementioned  criterion  (i.e.,  to  assume 
that  g/i0=f  and  determine  the  resulting  Kc  and  6  f),  we 
proceed  to  compare  the  scattered  field  from  surface  0D 
(b  =  2.27)  at  500  Hz,  with  predictions  from  the  RFM. 
Figure  1  illustrates  the  effect  on  the  RFM  and  H/K 
models  of  varying  the  slope  of  the  spectrum  (b)  in 
computing  S  f  from  Eq.  (5).  Note  first  that  for  this  case 


scattering  described  by  RFM  is  strongly  dependent  on 
the  spectrum  slope.  And  second,  the  best  fit  occurs  at 
about  b  =  2.14.  This  departure  from  the  expected  b  =  2.27 
is  not  particularly  troublesome  because  spectral  analysis 
of  the  actual  surface  realizations  resulted  in  spectral 
slope  values  between  2.10  and  2.35.  Examination  in 
greater  detail  bears  out  this  best  fit  value  for  b  and  also 
confirms  that  gMO  is  precisely  1.0. 

To  test  the  robustness  of  the  criterion  the  RFM  was 
run  for  surface  0D  at  three  other  frequencies  (250  Hz,  1 
kHz  and  2  kHz).  In  each  case  there  is  a  narrow  range  of 
b  values  that  fit  the  H/K  curve  within  3  dB.  The  value 
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Fig.  3.  Comparison  of  the  RKM  and  H/K  results  at  1  kHz  for  surface 
0D  assuming  that  j^0  =  1.0  and  b  =  2.14. 
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I'ig.  2.  Comparison  of  the  RK.M  and  H/K  results  at  250  11/  for 
surface  0D  assuming  that  =  1.0  and  b-2.14. 


I'ig.  4  Comparison  of  the  RKM  and  H/K  results  at  2fW 
surface  01)  assuming  that  =  1.0  and  b  -  2.14 
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TABLE  II.  BAND-LIMITED  RMS  ROUGHNESS  FOR  SURFACE 
OD  FOR  =1.0  CORRESPONDING  TO  THE  BEST  FITTING 
PARAMETERS  a  =  6.1(T3  AND  b  =  2.14. 


<V(m) 

<*f(deg) 

X 


250Hz 

5001b 

tkllz 

2kHz 

0.47 

0.23 

0.12 

0.06 

113 

13.1 

14.3 

15.3 

4.9 

5.4 

5.9 

6.4 

b=2.14  obtained  for  500  Hz  turned  out  to  be  a  good 
compromise  for  the  four  frequencies.  Figure  2,  3,  and  4 
show  the  results.  The  matching  between  the  RFM 
prediction  and  the  H/K  mode)  is  excellent.  Table  II 
gives  the  rms  height  of  the  microroughness  and  the  rms 
slope  angle  of  the  fine-scale  component  for  each 
frequency.  Also  shown  is  the  ratio  of  the  surface 
partition  wavenumber  to  the  acoustic  wavelength,  called 
the  partition  factor  x.  Notice  that  x  varies  signficantly 
and  consistently  over  the  three  octaves. 

The  same  procedure  of  looking  for  the  spectral 
parameter  b  that  leads  to  the  best  fit  for  gpo  =  1.0  was 
repeated  for  surface  other  surfaces  at  250  Hz.  Table  III 
shows  the  rms  values  of  the  fine-scale  slope  angle  and 
the  rms  height  of  the  microroughness  in  each  case.  In  all 
cases  the  criterion  gMO=  1.0  leads  to  the  best  fit.  The 
variations  seen  in  x  indicate  that  the  number  of  acoustic 
wavelengths  can  not  per-se  provide  a  criterion  to 
determine  the  partition  wavenumber. 


table  hl  best  fit  parameters  and  band-limited 
RMS  microroughness  for  surfaces  with  different 
degrees  of  roughness  at  a  constant  frequency 

OF  2S0  IIZ  (CUT-OFF  CRITERION:  =  1.0) 


surface 

label 

spectral  parameters 
a  b 

band-ltd  rms  values 
%  «»t  (<Rg) 

partttnn 

X 

OF 

2.3xl02 

1.75 

0.47 

10.3 

3.6 

01 

1.9xl02 

1.85 

0.47 

12.3 

3.5 

OB 

6.0xl03 

2.14 

0.47 

11.5 

4.9 

01 

4.8xl03 

2.18 

0.47 

10.6 

5.4 

OJ 

3.0xl03 

2.26 

0.47 

8.9 

7.5 

V.  CONCLUSIONS 

The  partition  wavenumber  between  the  two-scales  of 
roughness  in  composite-surface  scattering  theory  for 
surfaces  having  power-law  spectra,  typical  of  the  seafloor, 
has  been  investigated  through  numerical  simulation. 
Results  indicate  that  a  criterion  based  on  the  smoothness 
of  the  micro-scale  component  may  be  adequate  to 
establish  the  partition.  We  find  that  the  partition 
wavenumber  appears  to  be  set  by  the  physics  scattering 
to  be  fairly  precisely  g^o=1.0.  We  introduced  the 
hypothesis  that  a  facet  will  be  formed  such  that  its 
roughness  is  approaching  a  value  that  begins  to  destroy 
coherence  of  scatter  from  the  facet.  Larger  facets  are 
not  allowed  because  coherent  reflection  is  not  allowed. 
Smaller  facets  are  not  allowed  because  there  still  could 
be  a  coherent  reflection  from  a  large  one.  For  the  cases 
analyzed  here  the  hypothesis  appears  the  be  valid. 
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